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A new technique has been developed to estimate the probability that a nearby cloud to 
ground lightning stroke vvis within a specified radius of any point of interest 1 his process 
uses the bivariate Gaussian distribution of probability density provided by the current 
lightning location error ellipse for the most likely location of a lightning stroke and 
integrates it to determine the probability that the stroke is inside any specified radius of any 
location, even if that location is not centered on or even with the location error ellipse This 
technique is adapted from a method of calculating the probability of debris collision with 
spacecraft Such a technique is important in spaceport processing activities because it allows 
engineers to quantify the risk of induced current damage to critical electronics due to ncai by 
lightning strokes This technique was tested extensively and is now in use by space launch 
organizations at Kennedy Space Center and Cape Canaveral Air Force Station Future 
applications could include forensic meteorology 


Nomenclature 

0 = angle between the collision plane coordinates and the rotated collision plane coordinates 

p xz = correlation coefficient of x and z 
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o K = standard deviation of x coordinate of the diagonalized covariance ellipse in the rotated coordinate 

system 

<j/y = standaid deviation of z coordinate of the diagonalized covariance ellipse in the rotated coordinate 

system 

a x = standard deviation of x 

a z = standard deviation of z 

a x = standard deviation of x in the rotated coordinate system 

cr = standard deviation of z in the rotated coordinate system 

Vk = x coordinate of target circle in the (X Z ) coordinate system 

fi H = z coordinate of target cucle in the (X Z ) cooidinate system 

A = collision cioss sectional area (nautical miles 1 nmi 2 ) 

d0 = the angle between two points on the target ellipse 

dH = integration step 

H = intermediate variable in the lightning probability algorithm 

P = probability 

pdf = probability distribution function 

r A = radius of circle of cross sectional area A (nautical miles nmi) 

R = radius of ellipse (nmi) 

R1 = distance to first point on target ellipse (nmi) 

R2 = distance to second point on target ellipse (nmi) 

x = horizontal rectangular cooidinate in the collision plane (nautical miles mm) 

x = horizontal lectangular coordinate in the rotated collision plane (nmi) 

x = transformation variable that circularizes the x component of the probability ellipse (nmi) 

x c = nominal distance of closed approach of two colliding objects (nmi) 

W = intei mediate variable in the lightning probability algorithm 

= vertical rectangular coordinate in the collision plane (nautical miles nmi) 

= vertical rectangular coordinate in the collision plane (nautical miles nmi) 


z 



z 


transformation variable that circularizes the z component of the probability ellipse (nmi) 


I Introduction 

r J^Hb ability to accuiately estimate the probability that an individual neaiby cloud to giound lightning stroke was 
within a specified distance of any specified spaceport processing facility at Kennedy Space Center (KSC) or Cape 
Canaveral Air Force Station (CCATS) is important to processing payloads and space launch vehicles before launch 
Such estimates allow engineers to decide if inspection of electronics systems aboard satellite payloads space launch 
vehicles and ground support equipment is warranted due to induced currents from that stroke If induced current 
damage has occurred inspections of the electronics are critical to identify required Fixes and avoid degraded 
performance or failure of the satellite or space launch vehicle However inspections are costly both Financially and 
in terms of delayed processing for space launch activities As such it is important these inspections be avoided if 
not needed At KSC/CCAFS one of the main purposes of the Four Dimensional Lightning Surveillance System 
(4DLSS) [1 2] is detection of nearby strokes and determination of their peak current to support decisions to inspect 
electronics [3 5] The high frequency of lightning occurrence in East Central Tlorida combined with the large 
amount of complex sensitive electronics in satellite payloads space launch vehicles and associated facilities makes 
those decisions critically important to space launch processing The 4DLSS provides the data foi 50th peiccntilc 
location error ellipses for the best location for each stroke which is then scaled to 95th oi 99th percentile ellipses 
depending on customei requnements This error ellipse is necessarily centered on the best location of the lightning 
stroke The 4DLSS however has not been able to provide the probabil it> of the stioke being within a customer 
specified distance of a point of interest This paper presents a new method to convert the 4DLSS 50th percentile 
location enor ellipse for best location of any stroke into the probability that the stroke was within any radius of any 
facility at CCATS/KSC This technique could be adapted for use with National Lightning Detection Network 
(NLDN) data This new facility centric technique is a significant improvement over the stroke centric location error 
ellipses the 45th Weather Squadron (45WS) has provided in the past This technique is adapted from a method of 
calculating the probability of debris collision with spacecraft [6 8] 



II Methodology 


A Background 

In spacecraft collision probability and other applications at the instant of nominal closest approach the 
position uncertainty of the collision object relative to the asset being protected is described by a bivariate Gaussian 
probability density function (pdf) [6 8 9 10] as shown in the following equation 


/(*.*)= 



(i) 


where o x and o 2 = the standard deviations of x and z p xz = correlation coefficient of \ and / \ and 7 aic the 
designations for the rectangular coordinates in the collision plane 

The probability of collision (Lq 2) is given by the two dimensional integral where A is the collision cross 
sectional area which is a circle with radius r A [6] 


P = jjf(x,z)dxdz ( 2 ) 

A 

There is no known analytical solution to the above integral when the two standard deviations a x and o z 
are not equal The solution is found by performing a numerical integration of the two dimensional Gaussian pdf 
[6 8 10 ] 

The gcometiy used for spaceflight collision probability can also be used for estimation of the probability 
of an individual nearby lightning stroke contacting the surface within a specified distance of a specified point of 
interest as shown in Tig 1 Two methods of integrating the above probability weic tested and gave identical results 
The first solution method was based on an algorithm by Patera [8] and the second solution method was based on an 
algorithm by Chan [1 1] Chan s algorithm ran much faster and therefore was selected as the algorithm for the 4}WS 
lightning probability program 



B The numerical integration technique [11] 

This numerical integration technique is one in which the miss distance is given by a non central chi distribution 
with unequal variances [11] The covariance matrix corresponding to the bivanate Gaussian pdf in Eq 1 is [6] 


c = 


Px 


Px VxV 

a 1 


(3) 


When the correlation coefficient p xz is not zero there are undesirable off diagonal terms that ovei ly complicate the 
calculation In order to eliminate these terms the coordinate system (x z) is rotated to a new coordinate system 
(x z ) such that the majoi and minor axes of the ellipse associated with the covariance are aligned along the 
coordinate axes and the new covanance matrix is [6] 


C'= 


<*x 

0 


0 


(4) 


The angle 0 between the two coordinate systems is [6] 


0 = — tan' 
2 


2 P a o x c 




(5) 


The KSC/CCAhS 4DLSS system does not provide the covanance matrix but instead piovides the semi majoi 
axis semi minor axis and the orientation of the senn major axis of the 50% location enoi ellipse relative to north 
Therefoie the angle 0 in Eq 5 is found using geometry where 0 is the angle between the semi major axis of the 
lightning location unceitainty ellipse and line connecting the center of the lightning uncertainty ellipse and the 
center of the aica of interest 

In the (x z ) system the Eq I pdf becomes [6] 


/(*'»*)= 



2xo x o 


( 6 ) 


and Eq 2 the collision piobabihty becomes [6] 



where 


A'= A,r Al =r x'=x cos 0,z n -x sin# (8) 

’ A' A p e p e ' 

Tor spacecraft collision x c is the nominal distance of closest approach of the two colliding objects and (x p 7 p ) 

are the coordinates of the spacecraft relative to the debris Tor lightning strike probabi ! it\ \ c (the distance between 

the position of the centci of the strike location ellipse and the position of the target area) is calculated using the 

Haversine distance foimula 

The standard deviations in the new rotated coordinate system are calculated by dividing the semi major and 
semi minor axes of the 50% lightning positional confidence ellipse by the scaling constant used to scale standard 
error to the 50% confidence level The scaling constant is 

k = V-2* ln(l -0 50) (9) 

The probability is given by [1 1] 

P = —j=L — [L -{H-v, +e -(H + „ u fn a l L f(z ] ) + etf{z i )]dfj (13) 

2^2tco h l 

where 


V(»'-// , )+^]/V *7, 


(14) 


^2 = 


The parameters /u K and ju H are the coordinates of the target circle in the (X Z ) coordinate system and a K 
and <7 H are the standard deviations of the diagonalized covariance ellipse shown in Eq 4 The derivation of 


equations (13) and (14) above is shown in further detail in [1 1] A detailed example of the calculations using a real 
world case is provided in Appendix A 



in Evaluation 


The probability that any lightning strike is within any radius of any point of interest would be extremely difficult 
to estimate intuitively As a result given the high impact of the decisions on space launch operations the tool 
developed for this application was extensively tested Tests were conducted and are discussed in the following 
sections 1) known mathematical solutions and 2) examination of real uoild events The new technique passed all 
of the tests Tests were also conducted to assure the probabilities calculated using the algorithm of Chan [11] 
matched probabilities calculated using the algorithm of Patera [8] The probabilities calculated from the two 
algorithms were identical and thus are not shown here 


A Test Set 1 

The first set of testing compaicd the lightning sti ike piobability calculated using the 45 WS lightning sti ike 
spreadsheet (which uses an adaptation of the numencal integration algorithm by [11] to the conesponding circular 
probability fiom the CRC Handbook of Tables for Probability and Statistics [12] Table 1 shows the probability 
from the new numerical integration technique for various inputs and the corresponding correct probability from the 
CRC Handbook The values matched to within a tenth of a percent These errors in the final digit may be due to 
round off error 


Table 1 Calculated probabilit) vs CRC Handbook probability for various inputs 
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B Test Set 2 


The second type of testing analyzed six real world lightning strikes near Space Launch Complex 39A on 3 
August 2009 Figure 2 shows the spreadsheet used to generate the lightning report for those six strikes Additional 
data on three of these six stukcs aic m Table 3 T hese strikes were selected because the closest point on the lightning 
position unccitainty ellipse was within 0 45 nautical miles of Launch Complex 39A the key radius for assessing the 
need to inspect electronics for induced current damage to the Space Shuttle Figures 9 through 1 1 are Google Maps 
depictions of three of these six strokes The probabilities for a small area around a facility even for a nearby stroke 
may appear to be surprisingly low For example one strike just 0 65 nautical miles away (Figure 9) had only a 1 1% 
probability of being within the 0 45 nautical mile radius of Launch Complex 39A All calculated probabilities are 
consistent with these leal world events 

The KSC Clcctiomagnetic Environmental Effects (EEE) Panel requested six moic real world lightning strikes be 
investigated These wcic recently investigated lightning strikes near Launch Complexes 39A oi 39B where there 
was camera verification of the location of the strike The ECE Panel wanted to compare the icsults of the new 
facility centric probabilistic technique to these cases where the true answers were known unambiguously The data 
used for this analysis aie in Table 4 Both 4DLSS and National Lightning Data Nctwoik (NLDN) cases uctc 
examined depending upon which sensor system lecorded the stroke CGLSS strokes were obtained from 45 WS 
4DLSS The NLDN usually provided flash data so NLDN return stioke data weic pui chased as special StrikcNet 4 
repoits from Vaisala Coiporation This was done to match the return strokes routinely provided by 4DLSS Figures 
12 through 14 show the probability results from these cases As with the previous real world tests all calculated 
probabilities weie consistent with these additional real world events 


4 Vaisala Inc 2006 Vusah StrikcNet Information URL http //www vaisah com/filcs/StnkcNu Biochurt pdf 




Table 2 Input values used for scenarios shown in Figures 9 through 11 
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IV Other Applications 

The techniques and methods described in this paper clearly have application reaching far beyond the space 
program uses for which they were designed The list of potential applications is many and varied and would be of 
interest to anyone seeking information pertaining to probability of lightning strike locations such as the power 
industry aviation or any industry sensitive to electrical overloads This methodology is also applicable to forensic 
meteorology [13] wheie the question of whether lightning struck at or near a particular location is an issue in 
litigation This technique could be applied to NLDN or any other cloud to giound lightning detection system that 
provides location enor ellipses 


V Conclusion 

A technique has been developed to calculate the probability that any nearby cloud to ground lightning stiokc 
occurred within any ladius of any point of interest In piactice this provides the probability that a nearby lightning 
stroke was within a key distance of a facility rather than within the error ellipses centered on the stroke This 
process uses the bivariate Gaussian distribution of probability density provided by the current lightning location 
error ellipse for the most likely location of a lightning stroke and integrates it to determine the probability that the 
stroke is inside any specified radius This new facility centric technique which was tested extensively is much more 
useful to the space launch customers and has superseded the lightning error ellipse approach discussed in [3 4] 
Potential othei applications were briefly discussed in section IV 







Appendix 


This appendix is an example of calculating the probability of any lightning stroke with a known error ellipse 
being within a circle of any radius around any point It is provided to clarify the calculation process An example 
calculation is shown in Table 4 

This example is a leal world event from a lightning strike near the Space Shuttle launch pad 39A at 21 30 GMT 
on 14 Oct 2009 (ref 1 igure 12) Although the lightning data usually are from the cloud to ground component of the 
Four Dimensional Lightning Surveillance System (CG 4DLSS) (Murphy et al 2008 Roeder 2010) in this example a 
lightning stroke from the NLDN is used We sometimes use StrikeNet reports that provide stroke data from the 
NLDN to double check the CG 4DLSS report 

1) Location of Launch Pad 39A 28 60827486 N (or 0 499309 radians) 

80 6041 1653 (i c 80 6041 1653 W or 1 406807 ladians) 1 his is also the ccntei of the circle in which the 
lightning piobability will be calculated 

2) Desned Radius Tor Probability of Lightning Aiound 39A 0 45 nautical miles 

3) Lightning Stioke Data Time/date 0235 GMT 16 August 2009 latitude 28 6069 N (or 0 499285 
radians) longitude 80 6087 (le 80 6087 W or 1 406887 radians) polarity/peak cuirent 48 0 kA semi 
major axis of 50% confidence location ellipse 0 6 km semi minor axis of 50% confidence location 
ellipse 0 4km orientation angle of location ellipse 82 (clockwise from north) 

Table 3 Lightning strike probability calculation process 

Step Action Equation and other information Calculation and Result 


1 Convert semi 1 nmi = 1 852 km 
major and semi 
minoi axes from 
km to nmi 


0 6 km = 0 324 nmi semi 
major axis 

0 4 km = 0 216 nmi semi 
minor axis 


2 Calculate 

distance from 
lightning stroke 
to centei of 
circle 


Haversine Distance formula 
Distance = Earth Radius * C 

• Earth Radius = 3443 920086 nmi 

• C = 2* Atn2[y/] — A -Ja ) 

Atn2 is a two parameter arc tangent function which 
leturns values in all four quadrants 

• A = sin(dlat/2)*sin(dlat/2) + cos( (target 
lat))* cos((stroke lat)) 
*sin(dlon/2)*sin(dlon/2) 

• dlat = latitude difference from target to 
stroke = 


Distance = 3443 920086 * * 
7 42 1 7\ 10 5 
= 0 2556 nmi 

C = 2*Atn2{sqi(l 
1 377x1 0 9 ) 
sqr(l j77\10 9 )} 

= 7 42 1 7\ 1 0 5 

A = sin(1 200x10 5 /2)* 
sin(l 200x1 0 5 /2)* 

+ 





3 Perform 
coordinate 
system rotation 
to eliminate the 
off diagonal 
term in the 
covariance 
matrix 


4 Calculate the 
standard 
deviations in the 
new rotated 
coordinate 
system 


28 60827° 28 6069° = 2 39959x10 5 
(radians) 

• dlon = longitude difference from circle to 
stroke = 80 60411° 80 6085° = 

7 99967X 10 5 (radians) 


cos(0 4993)*cos(0 4993)* 
sin(4 000x1 0 5 /2) 

*sin(4 000 x1 0 5 /2) 

= 1 3770x10 9 


• X = (Longitude of Target - Longitude of Stroke 

) * Cos (Latitude of Strike) 

• Z = Latitude of Target - Latitude of Stroke 

• 0 = a - ((tt/2) - Atn2(X Z) 

• a is the orientation angle of the 50% lightning 

positional confidence ellipse 

• X = miss distance*Cos(0 (coordinate system 

rotation angle)) 

• Z = miss distance*Sin(0) 


X = ( 1 4068- ( 1 4069) * 
Cos (0 4993) 

= 7 0231 X 10 5 
Z = 0 49931 -0 49929 
= 2 400 X I0 5 

0 = I 431 - I 571 - 
Atn2(7 023 X 10 5 2 400 
X 10 5 ) = 0 1896 

X = 0 2556*Cos(0 1896) 
= 0 2510 

Z = 0 2556*Sm(0 1896) 

= 0 0482 


• o x = Semi major axis of the 50% lightning 

positional confidence ellipse /elliptical scaling 
constant used to scale standard error to the 50% 
confidence level 

• o z = Semi minor axis of the 50% lightning 

positional confidence ellipse /elliptical scaling 
constant used to scale standard error to the 50% 
confidence level 

• Elliptical scaling constant k is 

yj-2* !n(l - probability) 


k = ij-2 * ln(l - 0 50) 

= 1 1774 

o x = 0 3240/11774 = 

0 2752 

g 2 = 0 2160/1 1774 = 

0 1834 
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Calculate the 
probability that 
lightning stroke 
was within the 
target aiea of 
intei est b> 
pei foiming a 
numerical 
integration using 
Simpson s rule 
of the lightning 
uncertainty 
ellipse over the 
area of the circle 
around the target 
of interest 


p = 






/[• 


(« M Y / 2a 


+ e (" / ) /2 ° j[tr/(Z,)+ei/(7 2 )]c/// 


Z| = 
= 


#-// 2 )-/'x]/V la K 

Vfr-// 2 ) + ^]/V2a A 


• p k and pn are the coordinates of the target circle in 
the (X Z ) cooidinate system 

• o k and o H are equal to o x and o z which are the 
standard deviations of the diagonalized covariance 
matrix 

• W = Radius around target 2 

• N = the number of iterations to perform in the 

integration (for this example N is set to 200) 

• DH - yfw / N = integration step 

• B C and D are intermediate variables in the 
algorithm corresponding to various parts of the 
probability equation shown above 


• B = V2 * o x 


• C = X /B 


W = Radius around target 7 
W = 0 45 2 = 0 2025 

DH = yfw / N 

Dll = V0 2025 /200 = 
0 00223 

B = V 2 * o x 

B= VI * 0 2752 = 0 3891 
C = X / B 

C = 0 2510/0 3891 = 

0 6451 

D = 1 / (2 * -Jin * o z ) 

D = 1/(2 * V2/r * 0 1834) = 

1 0874 


1 1 = iteration no * DH 



• A H zl z2 E r Erf(zl) Erf(z2) Q and sum arc 

intermediate variables in the algonthm 
corresponding to various parts of the probability 
equation shown above 

• A loop is performed forj = 1 to 199 This 

example is shown for j = 199 

• Sum = 0 

• Begin Loop here H =j* DH 

• A = 

• zl =A/B-C 

• z2 = A/B + C 

2 v 2 

• Erf(x) = error function = —j= \e~‘ dt 

v it o 

• E = (H - Z ) 1 / (2 * o z ’) 

• F = (H + Z ) 2 / (2 * o z 2 ) 

• QG) = (e E + e F ) * (brf(Zl) + Erf(Z2)) 

• sum = sum + (3 ( l) J )*Q0) 

• sum = sum + Q(0) + Q(N) 

• Probability = D*sum*DH/3 


QO) = (e E + e F ) * 
(Erf(Zl) + Erf(Z2)) 

Q( 1 99) = (e 2 372 + e 3 65 *)* 

( 0 5461+0 7179) = 

2 0467 X 10 2 

sum = sum + (3 ( 1) J ) * 

QO) 

sum = sum + (3 (1) '") 

* 2 047 X 10 2 
sum = 844 8952 
Lnd Loop 

sum = sum + Q(0) + Q(N) 
sum = 844 8952 + 2 9361 
+ 0 = 847 8317 

PiobabilUy = D * sum * 
DH / 3 

Probability = 1 0874 * 
847 8317 * 0 00225 / j = 
0 6914 




A = 4 494 X 10 ’ 


zl = A/B-C 

zl = 4 494 X 10 2 / 0 3891 

-0 6451 

zl = 0 5296 

z2 = A / B + C 

z2 = 4 494 X 10 2 / 0 3891 

+ 0 6451 

z2 = 0 7606 

Erf(Zl) = 
ErrorFunction(zl) 

Ei f( 0 5296)= 0 5461 

Erf(Z2) = 

Crroi 1 unction(/2) 

Eif(0 7606) = 0 7179 
E = (H - Z ) ’ / (2 * o z 2 ) 

E = (0 4478- 
0 0482) , /(2*0 1 834’) 

E = 2 372 

r = (H + Z) 1 /(2*o z ’) 
r = (0 4478 + 

0 0482) 2 /(2*0 I834 7 ) 

T = 3 654 



E = (H - Z ) 2 / (2 * c z 2 ) 

E = (0 4478 - 

0 0482) 2 /(2*0 1834 2 ) 

E = 2 372 

T = (H + Z) 2 /(2 *o z ’) 

1 =(0 4478 + 

0 0482)’/(2*0 1834’) 

T = 3 654 

Q0) = (c f + c ') * 
(Crf(Zl) Lrf(Z2)) 

Q( 1 99) = (e 2 372 He 1651 )* 
( 0 5461+0 7179) = 

2 0467 X 10 2 

sum = sum + (3 ( 1) J ) * 

QO) 

sum = sum + (3 (1) '") 

* 2 047X10 2 
sum = 844 8952 
End Loop 

sum = sum + Q(0) + Q(N) 
sum = 844 8952 + 2 9361 
+ 0 = 847 8317 

Probability = D * sum * 
DH/3 

Probability = 1 0874 * 
847 8317 * 0 00225 /3 = 
0 6914 
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Fig 1 Schematic diagiam of the angles used in probability calculation for a sample lightning location error 
ellipse a is the heading of the semi major axis of the lightning location uncertainty ellipse from true north 0 
is the angle between the semi major axis of the lightning location uncertainty ellipse and line connecting the 
center of the lightning uncertainty ellipse and the center of the area of interest 
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Fig. 2 Sample of lightning strikes where the closest point on the lightning position uncertainty ellipse was 
within 0.45 nmi of Launch Complex 39A on 3 August 2009. 




Fig. 3 Google Maps visualization of the 99% confidence uncertainty ellipse for one of the closest lightning 
strikes to Complex 39A on 03 August 2009. The black ellipse is the 99% lightning location uncertainty 
ellipse. The white circle is a 0.45 nmi radius around the point of interest. The center of the ellipse was within 
the 0.45 nmi radius. There is a 53.8% probability that the lightning occurred within that radius. 
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Fig. 5 Google Maps visualization of the 99% confidence uncertainty ellipse for nearby lightning strike to 
Complex 39A on 03 August 2009. The black ellipse is the 99% lightning location uncertainty ellipse. The 
white circle is a 0.45 nmi radius around the point of interest. Figure 11 shows a probability of 1.1% of the 
lightning strike occurring within the 0.45 nmi radius. 



Fig. 6 Illustrates a probability of 69.1% of a lightning strike of amplitude -43.0 kA detected by NLDN 
occurring 0.26 nmi from the center of Launch Complex 39A on 8/16/2009. The black ellipse is the 99% 
lightning location uncertainty ellipse. The white circle is a 0.45 nmi radius around the point of interest. 
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